Thin topological insulator film in a perpendicular magnetic field 
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We report on a study of an ultrathin topological insulator film with hybridization between the 
top and bottom surfaces, placed in a quantizing perpendicular magnetic field. We calculate the full 
Landau level spectrum of the film as a function of the applied magnetic field and the magnitude of 
the hybridization matrix element, taking into account both the orbital and the Zeeman spin splitting 
effects of the field. For an undoped film, we find a quantum phase transition between a state with a 
zero Hall conductivity and a state with a quantized Hall conductivity equal to e 2 /h, as a function of 
the magnitude of the applied field. The transition is driven by the competition between the Zeeman 
and the hybridization energies. 
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I. INTRODUCTION 

Topological insulator (TI) is a new phase of matter, 
which has recently been discovered in materials, long 
known and extensively studied in the fields of traditional 
semiconductor physics (HgTe) and the physics of thermo- 
electrics (Bi 2 Se3 and Bi 2 Te3)ji~— Somewhat uncharacter- 
istically of the "standard" order of things in condensed 
matter physics, TI phase was first predicted to occur in 
these materials theoretically^ and only later seen exper- 
imentally— One obvious reason for this is that the main 
experimental manifestation of the TI phase, namely the 
occurence of metallic edge states in insulating samples, 
is a somewhat subtle effect, unlike, e.g., the spectacu- 
lar quantization of the Hall conductivity in a close rela- 
tive of TI, the two-dimensional electron gas (2DEG) in 
the quantum Hall effect (QHE) regime. While theoret- 
ically TI phase can in fact be characterized by a quan- 
tized physical quantity, the topological magnetoelectric 
susceptibility;^ measuring this quantity experimentally 
is not straightforward. Currently, the only realistic pro- 
posal involves measurement of the magneto-optical re- 
sponse in a thin TI film when the time-reversal symme- 
try is broken, with the Faraday and Kerr rotation angles 
predicted to acquire universal values in this case£~— 

As far as transport phenomena in TI are concerned, 
theoretical effort has mostly been concentrated on study- 
ing the effects of the characteristic spin-momentum lock- 
ing of the helical TI surface states, which has an obvious 
potential in spintronics. 9-18 An important practical is- 
sue, relevant to all transport studies of TI, is that all cur- 
rently known TI materials, while theoretically insulators, 
in reality are metallic in the bulk. This happens due to 
unavoidable doping, introduced by impurities and crys- 
tal lattice defects. One possible way to deal with prob- 
lem is to grow TI samples in the form of thin filmsjiSr— 
This allows to both directly reduce the bulk contribution 
to transport by simply reducing the bulk volume and to 
control the carrier concentration in the bulk through gat- 
ing iS^r. 2 ^ 

Ultrathin TI films also have interesting physical prop- 
erties, such as the already mentioned above universal 



magneto-optical response^— possible excitonic superflu- 
idity^ strongly improved thermoelectric performance, 27 
quantum spin Hall 28-29 and quantum anomalous Hall 
(QAH) effects,— which have so far been explored only 
theoretically. Of particular interest are the effects which 
arise in the ultrathin limit, when the top and bottom sur- 
faces of a TI film start to hybridize. •■— We have recently 
pointed oul^ that magnetic response, in particular re- 
sponse to field, applied in the plane of the film, is highly 
nontrivial in this limit, with the film undergoing a topo- 
logical insulator to semimetal quantum phase transition 
function of the field. 
In this paper we study the properties of TI thin films 
in a perpendicular quantizing magnetic field. We calcu- 
late the full Landau level (LL) spectrum of a TI film with 
hybridization between the top and bottom surfaces and 
the Hall conductivity. One would normally expect the 
Hall conductivity to be zero in an undoped charge-neutral 
film, which is indeed what happens when hybridization 
is stronger than the Zeeman splitting. As the magnetic 
field is increased, however, we find a quantum phase tran- 
sition, at which the Hall conductivity jumps from zero 
to e 2 /h. The transition happens when the Zeeman en- 
ergy becomes larger than the hybridization energy and is 
shown here to be a direct consequence of the well-known 
characteristic feature of the Dirac-fermion LL spectrum, 
namely the existence a zero-energy LL. The transition 
we find is analogous to the transition to the QAH state 
discussed in Ref.[3(| in the context of TI films, doped 
with magnetic impurities. Here we demonstrate that the 
same physics can be realized in a more straightforward 
and currently experimentally accessible way, simply uti- 
lizing an applied magnetic field. Indeed, experimental 
work on magnetotransport properties of thin TI films in 
quantizing perpendicular magnetic field has already ap- 
peared in the literature^r— 



II. LANDAU LEVEL SPECTRUM OF AN 
ULTRATHIN TI FILM 

We start from the following Hamiltonian of an ultra- 
thin TI film in the presence of a perpendicular magnetic 
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field: 

H = J drtf^r) v f t z (z X er) • (-iV + ^a) + A z a z 

+ A t T x ]V(r). (1) 

Here vf is the surface Dirac cone Fermi velocity, A z = 
g/is-B/2 is the Zeeman energy, associated with the ap- 
plied magnetic field B = Bz (we will assume B > 
henceforth), A t is the hybridization matrix element (we 
will also assume A t > 0) and Ti = 1 units are used. We 
have introduced Pauli matrices er and r to describe the 
real spin and the which surface pseudospin degrees of 
freedom and suppressed the explicit spin and pseudospin 
indices for clarity of notation. 

To diagonalize Eq.(p}, we choose Landau gauge for the 
vector potential A = xBy and define LL ladder operators 
in the standard way as: 

a = {£/V2)(tc x - nr y ), ^ = {£/V2)(tt x + iv y ), (2) 



where 



where I — ^Jc/eB is the magnetic length and tv = 
— iV + (e/c)A is the kinetic momentum. The single- 
particle Hamiltonian operator can then be written in 
terms of the ladder operators as: 



H 



■■{cr+a- 



a a 



A z a z + A t r* 



(3) 



where we have introduced a characteristic frequency 
lob = vp/£, which plays a role, analogous to the cyclotron 
frequency in the LL spectrum of a regular 2DEG. 

It is clear from Eq.((3]) that the single-particle eigen- 
states have the following general form: 



\nas) = u 
+ u 



as I 
nTfl 
as 



l,T,t> + <^|n,T,i) 
l,B,t)+< s Bi \n,B,l). 



(4) 



Here |n,T(B),t (J,)) is the n-th LL eigenstate on the 
top (bottom) surface with spin up (down), a = 0, 1 and 
s = ± label the four eigenstates of Eq. © , corresponding 
to each LL index n — 0, . . . , oo, and u" s are the corre- 
sponding complex four-component spinor wavefunctions. 
We have also suppressed the intra-LL orbital label for 
brevity. The physical meaning of the a and s indices will 
become clear shortly. 

The problem of finding the eigenstates of Eq.Q thus 
reduces to diagonalizing the following 4x4 matrix: 



/ A, 

V o 



-A z 


A, 



A t ^ 

A t 

A z iuisV^n 



(5) 



Diagonalizing ([5]), we find the following LL spectrum: 



e na ± = (-lT^2uln + (A, ± A t f. (6) 
The corresponding spinor wavefunctions are given by: 

~ 1) /nas+) ^/nas-j^( 1) fnas+ifnas—]: 

(7) 



fnas± — ~ J 1 ± — ^ 
Z V £r, 



(8) 



The LL spectrum thus consists of two sets: electron-like 
(e„o±) and hole-like (e„i±) Landau levels. Within each 
set, every LL is further split into a doublet (±). Note 
that this splitting results from the presence of both the 
Zeeman spin splitting and the hybridization, and van- 
ishes if any one of them is zero (except for the n = LL, 
which is split whenever A t is nonzero). 

The n = LL are special. As seen from Eq. ©, elec- 
trons in these levels are fully spin polarized, i.e. only 
the spin-down states are occupied. Correspondingly, the 
n = level is only split into two sublevels, unlike all 
other LL, which are split into four. This is the well- 
known "zero-mode anomaly" of the Dirac-fermion LL 
spectrum;^— which will play an important role in our 
story. 

Explicitly, the two n — LL wavefunctions are given 
by: 
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1 



= —(0,1,0,1), «J+ = -=(0,-1,0,1), (9) 



V2 



V2 V 



when A z < A t . The corresponding LL energies in this 
case are: 

e 00 _ = \A Z - A t ], eoi+ = -|A, + A t |, (10) 
On the other hand, when A z > At, Uq~ is replaced by 



ul~ = -=(0,1,0,1), 
with the corresponding energy given by: 



eoi- 



|A, - A t 



(11) 



(12) 



Thus one of the two n = sublevels moves from the 
set of electron-like levels into the set of hole-like levels. 
This change of the character of the LL spectrum as a 
function of A z is illustrated in Fig.Q] As we demonstrate 
below, the transition to purely hole-like n = LL is 
manifested in a jump of the Hall conductivity at fixed 
chemical potential ep = by e 2 /h. 



III. HALL CONDUCTIVITY 

Given the above LL spectrum, we can now evaluate 
the Hall conductivity of the TI film. Assuming the Fermi 
level to be always in a gap between the LL, the Kubo- 
formula expression for the Hall conductivity reads: 



? 2 

2vr 4^ 

nas^n ol s' 



x (n'a's'\T z a x \nas)] 



Im {(nas\T z cr v \n'a's') 

nF(inas) — np^n'a' 



(£,. 



(13) 
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ble for all n > LL, this term gives: 



a) 



w VJ\L 




b) 



FIG. 1. (Color online) LL density of states of a thin TI film 
(broadening added by hand). Dashed line represents the zero 
of energy, set at the Dirac point position in the absence of 
the magnetic field, (a) A z = 0, At > case. Only the n = 
LL is split and the spectrum is particle-hole symmetric, (b) 
A z > At > case. Both n — LL are hole-like and particle- 
hole symmetry of the spectrum is broken. All LL splittings 
in the figure are exaggerated relative to the main LL spacing. 



where Hp is the Fermi-Dirac distribution (we will as- 
sume zero temperature in all calculations). The matrix 
elements in Eq. (jl3p are nonzero only when s — s' and 
n! = n ± 1. Evaluating the matrix elements, we obtain 
the following expression for the Hall conductivity: 



•> — X! X! ^ {( 2n + !) [ n F^nas) ~ n F (e n+las )} 

n=0 ct=0,l s=± 

n F (e ) n F (e nas ) 



4tt 

(A, + sA t ) 



(14) 



>:ry 



— (2/1 + 1), n > 1, 



(15) 



where we have restored h for clarity. This result is the 
same (up to a factor of 1/2, accounting for half the degen- 
eracy) as the well-known result for the Hall conductivity 
of graphene^ 

The contribution of the n — LL is contained in the 
second term in Ea. (fl4l . which also contains the physics, 
unique to TI thin films. The effect of this term is most 
easily seen at the point e F = 0, which corresponds to 
the charge-neutrality point in the limit A z = (A 2 > 
breaks particle- hole symmetry, as seen from Fig. [1}. In 
this case we obtain the following expression for the Hall 
conductivity (note again that A z , A t > 0): 



2/? 



[sign(A 2 - A t ) + 1] 



(16) 



Thus one can see that cr xy , evaluated at fixed Fermi 
energy e F — 0, jumps from (A z < A t ) to ±e 2 /h 
(A 2 > A t ) as a function of the applied magnetic field. 
This transition is the result of the change of the charac- 
ter of the n = doublet of LL, which happens as a func- 
tion of A z /At, and which was described in Section ITll 
To reiterate, when A z < A t , one of the n = LL (u°~) 
is electron-like, while the other one (u 1+ ) is hole- like. 
When A 2 > A t , both n = LL become hole-like (w 1± ), 
and the Hall conductivity thus increases by e 2 /h. This 
change in the value of the Hall conductivity is of course 
only observed when the Fermi level is initially within the 
hybridization gap of the n — LL. 

It is important to note that the above transition at 
a fixed Fermi energy is associated with a small in- 
crease of the (two-dimensional) electron density by 5n e ~ 
l/27r£ 2 |A,=A t = mAt/n, since the high-field side of the 
transition A z > A t differs from the low-field side by an 
extra filled LL. If the charge density is held fixed and 
one starts from the low-field limit when the Fermi level 
is in the hybridization gap between the two n = Lan- 
dau sublevels, no transition occurs upon increasing the 
field, as the Fermi level moves with the field and always 
stays in the gap. However, the fixed Fermi energy sit- 
uation might in fact be more physically relevant for a 
sample in a transport-measurement setup, attached to 
external conducting leads (admittedly, this is a subtle is- 
sue and what actually happens may depend on details of 
the experimental setup and various characteristics of the 
sample) . 



IV. DISCUSSION AND CONCLUSIONS 



The first term in Eq. (jl4l corresponds to the well-known 
expression for the Dirac-fermion Hall conductivity. As- 
suming n electron-like LL are filled and neglecting the 
Zeeman and hybridization splittings, which are negligi- 



In this paper we have studied the properties of an ul- 
trathin TI film, in which the top and bottom surfaces are 
hybridized, placed in a perpendicular magnetic field. We 
have calculated the LL spectrum of such a film, taking 
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into account the Zeeman spin splitting. We have also 
evaluated the Hall conductivity of the film and demon- 
strated that the Hall conductivity has a nontrivial de- 
pendence on the ratio of two energy scales: the Zeeman 
energy A z and the hybridization energy A t . Namely, we 
have shown that a thin film with hybridization between 
the top and bottom surfaces undergoes a quantum phase 
transition from a "trivial insulator" state, in which the 
Hall conductivity at charge neutrality is zero to a "Hall 
insulator" state, in which the Hall conductivity is equal 
to e 2 /h. The transition happens when A 2 — A t and is 
associated with a nonanalytic contribution to the Hall 
conductivity, which comes entirely from the n = LL. 
The transition is accompanied by a change in the charac- 
ter of the n — LL. On the low-field side of the transition 
one of the n = sublevels is electron-like, while the sec- 
ond one is hole-like. On the high-field side both LL are 
hole- like. 

Finally, let us briefly discuss experimental observabil- 
ity of the proposed effect. To this end we estimate the 
energy scales, associated with the TI film LL. Taking 
B ~ 10 T, we obtain A z ~ 1 meV (using g ~ 1 for a 
very conservative estimate), and uib ~ 10 meV. Assum- 
ing A t ~ A z (which is needed to observe the phase tran- 
sition we have described), we then have > A t ,A z . 
Thus observing the "trivial insulator" to "Hall insula- 



tor" phase transition might require more stringent con- 
ditions than a general observation of the QHE in this sys- 
tem would require. One needs low, i.e. less than about 
10 K temperatures, and a clean enough sample, that a 
~ 1 meV hybridization gap between the two n = LL is 
not washed out by disorder. However, recent experiments 
actually suggest a much larger value of the ^-factor for 
the surface states in Bi2Se3 and related materials ; 32 ^ 3 
as large as g w 50 (which is not unexpected due to the 
strong spin-orbit interactions in the bulk material^), and 
this makes the proposed phase transition much more eas- 
ily observable. But even the most stringent conditions, 
required if g ~ 1, are in fact standard in semiconduc- 
tor 2DEG QHE measurements and thus should be easily 
achievable, at least for Molecular Beam Epitaxy grown 
TI films. 
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